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<^ • Abstract 

Starting from exact expression for the dynamical spin susceptibility in the time-dependent den- 
sity functional theory a controversial issue about exchange interaction parameters and spin-wave 
excitation spectra of itinerant electron ferromagnets is reconsidered. It is shown that the origi- 
nal expressions for exchange integrals based on the magnetic force theorem (J. Phys. F 14 L125 
(1984)) are optimal for the calculations of the magnon spectrum whereas static response function 
is better described by the "renormalized" magnetic force theorem by P. Bruno (Phys. Rev. Lett. 
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■^ • 90 087205 (2003)). This conclusion is confirmed by the ab initio calculations for Fe and Ni 

o' 

t^ ■ PACS numbers: 75.30.Et 75.30.Ds 71.15.Mb 75.50.Bb 

a 

o 
o 

> 

X 



An efficient scheme for the first-principle calculations of exchange interaction parame- 



ters in magnets based on a so-called "magnetic force theorem" (MFT) 



m 



in the density 



functional theory is frequently used for analysis of exchange parameters for different classes 
of magnetic materials such as dilute magnetic semiconductors J3|, molecular magnets j^, 
colossal magnetoresistance perovskites [5J, transition metal alloys |6|, hard magnetic ma- 
terials such as PtCo |7[ and many others. Recently this method was generalized to take 



into account the correlation effects and_successfully used for the quantitative estimation 
of exchange interactions in Fe and Ni 



|9|. At the same time the formal status of this 
approach is still not well-defined since a general mapping of formally rigorous spin density 
functional to an effective classical Heisenberg Hamiltonian can be done only approximately. 
It was noticed already in the first work on the MFT |jj that only the expression for spin 
wave stiffness constant D is reliable. In terms of the diagrammatic many-body approach 
it means that the exchange integrals (J,,) in general should contain vertex corrections ICrl 

fl Q R 

which are neglected in the simple expression [Jj, |2] (see Ref. |2l). At the same time, using a 
general expression for spin- wave stiffness due to Hertz and Edwards 11| one can prove that 
the vertex corrections to D are cancelled for any local approximation for the self-energy (or, 
in the density functional method, for the local exchange correlation potential) [ij]. Recently 
P. Bruno has suggested [13| corrections to the MFT and consequently to the expressions for 



Jij (see also Ref. |1J). It is important to note that, first, the new expression for D coincides 
with the old one and, second, that these corrections for the case of itinerant electron mag- 
nets are formally small in adiabatic parameter 77 = w/A where ccj is a characteristic magnon 
frequency and A is the Stoner spin splitting. At the same time, mapping of the local spin 
density approximation (LSDA) onto the classical Heisenberg model itself is valid only in the 
adiabatic approximation r/ — > [ij, |2, ll6| . If we are interested in higher order effects in the r] 
it might be needed different effective exchange parameters for different physical properties. 
We will show on few examples that this is exactly the case. It turns out that the spin-wave 
excitation spectrum should be calculated in terms of "old" exchange integrals 111 M whereas 
for static properties "new" exchange integrals [IJ, H^ are more appropriate. 

The most reliable way to consider spin-wave properties of itinerant electron magnets 
in the framework of the spin density functional theory is the use of frequency dependent 
magnetic susceptibility [l7|, |l8|, |l9|. One should start from the time-dependent density 



functional theory in the adiabatic approximation (ADA-TDDFT) J2y, |2l|. We proceed with 



the Schrodinger-like equation within the self-consistent ADA-TDDFT potential 

H = -V^ + V{v) - i(B,,(r) + B,,,(r))a (1) 

(Slater's units are used here) where V^(r) is an effective potential, 'Bextij) and 'Q^cij) are 
external magnetic field acting on spin and exchange- correlation field, respectively. The 
adiabatic approximation means that the functional dependences of exchange-correlation po- 
tential and field on the charge and spin density are supposed to be the same as in the 
stationary case. In LSDA one has 



J r — r on 



J |i — 1 I on 
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B,, = -2—j—[ne,,] (2) 
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where n and m are charge and spin density, Exc is the exchange-correlation energy per 
particle, Vext is the external potential, i.e. the Coulomb potential of nuclei. To calculate 
the spin susceptibility we will assume ^exti^^) ^ 0. It leads to the effective complete "non- 
equilibrium" field 

SB^t = SB:,t + ^^^^ (3) 

where a(3 are Cartesian indices and the sum over repeated indices is assumed. 

By definition of the exact non-local frequency-dependent spin susceptibility x"^ the vari- 
ation of the spin density is equal to 

6m- = r^6Bi, (4) 

Operator product is defined here as usual: 

{X^){r) = j dr'x{r,r'Mr') (5) 



201 that in the time-dependent density functional theory 



On other hand, it was shown in Ref. 
we should have exactly 

<5m- = To'SBt, (6) 

where Xo is the susceptibility of an auxiliary system of free Kohn-Sham particles. It has 
been demonstrated originally for the charge excitations but the generalization of TDDFT 



for the spin-polarized case |21|, shows that the equation (jH)) holds also for spin excitations. 
Comparing two expressions for 5m" we have the following equation 

r' = xf + Xo'^-^x'' (7) 

This "RPA-like" equation is formally exact in ADA-TDDFT. For the local spin density 
approximation (Eq.Q) one has 

5m^ m \ ^ m? ) m? dm 

The first term in Eq.(jHJ is purely transverse and the second one is purely longitudinal with 
respect to the local magnetization vector. In collinear magnetic structures there are no 
coupling between the longitudinal and transverse components and for the transverse spin 
susceptibility we have the following equation: 

X+-(r,r',a;) = Xo+-(r,r',.;) + y"drV(r,r",a;)/..e(r")x+-(r",r',.;) (9) 



where 
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Bxc ^^Q^ 



is an exchange-correlation "Hund's rule" interaction, the magnetic and charge electron den- 
sity being defined as usual 



m 
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The bare susceptibility has the following form: 



Xo+-(r,r',a;) = ^ ^'' f^^ r,,{r)ij^i{r)r.d^')i^,l{r') (12) 



where ■j/'^o- and e^o- are eigenstates and eigenvalues for the Kohn-Sham quasiparticles 
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Ho = -V + V{r) (13) 



Ho - -crB-,^ ] i>^^ = e^^i)^^ 



and f^a = f (e^a) is the Fermi distribution function. 

Although the longitudinal spin susceptibility is not necessary to consider the exchange 
interactions it is instructive to write also an explicit expression for it. The derivation is 



similar to those presented above with a small complication, since we have to consider sep- 
arately the response of spin-up and spin-down electrons. Suppose we have an external 
perturbation dV^^f-. It leads to the change of the exchange correlation potential V^^ = g^ 

{n„ = i (n + am)), namely, 

5V:, = U,„>5n,., (14) 

^aa' — t: t: 
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which gives the total perturbation ^V^^t = ^^ext + ^^xc- O^^ can introduce the response 
functions i^°"°" by the formal expression 

dn'^ = K^^'SV::,, (15) 

Then the longitudinal spin susceptibility can be expressed in terms of the i^-functions as 

X'^ = 1 (i^TT + ;^U _ j^n _ ;^iT) (16) 

At the same time, according to the general scheme of TDDFT, we have 

5n^ = X,5VZt (17) 

where 

^.(r,r') = E;7^7^T^CW^..W^..(r')C(r') (18) 

Comparing these two expressions for 5n'^ one obtains 

K^^ = X^+X^U^^K^^+X^U^^K^^ 

/sfti = X^U^iK^^ + X^U^^K^^ 

K^^ = X^U^^K^^+X^U^K^^ (19) 



Similar expressions have been obtained in the RPA for the Hubbard model in Ref |22. A 
coupling between the longitudinal spin and density degrees of freedom is important also for 
the electronic structure calculations which take into account correlation effects (fll l23l| . 

Let us continue the derivation of useful expression for the transverse susceptibility 
(Eq.(jHl)). In order to consider the case of small u it is useful to make some identical transfor- 



mations of the kernel (|T^ similar to the Hubbard model consideration 2^. Using Eq.()13|) 



one can find 

B^^ifj^^rui = i^-i - ^a«t) Ci^MT + V(V^;.T VCi - CiVV^^t) (20) 

Substituting Eq.(I2ni) into Eq.dHI) we obtain 

{xt'B^c){r, r', uj) = m{r)S{r - r') - ujxt'i^, r', uj) (21) 

where we used the completeness condition 

J2r,Ar)%Ar') = S{r-r') (22) 

Substituting Eq. (j2H) into Eq. (jl2j) we can transform the latter expression to the following 
form ^ 

X =Xo +Xo )C =X6 +X - ^Xo —X + —X^ (23) 

m mm 

or, equivalent ly, 

X+- = m[uj-{xo^-y'Ay' (24) 

where 

A(r,r',.;) = E ^T !^t r,,ir)Ui^)'^ [V^MT(r')VCi(r') " Ci(r')V^,T(r')] (25) 

Using Eqs.(IT21),(I21I) one has finally 

X+-= (m + A) (uj-I,,aY' (26) 



which is exactly equivalent to Eq. but much more suitable for investigation of the magnon 
spectrum. Spin wave excitations can be separated from the Stoner continuum (e.g., para- 
magnons) only in the adiabatic approximation, which means the replacement A(r, r',u;) by 
A(r, r', 0) in Eq. ()26|) . Otherwise one should just find the poles of the total susceptibility, and 
the whole concept of "exchange interactions" is not uniquely defined. Nevertheless, formally 
we can introduce the effective exchange interactions via the quantities 

l](r,r',cu)=4eA(r,r',cu). (27) 

Substituting Eq.(I2ni) into Eq.(j2ni) we get 

A(r, r', uj) = J2 ^^rf^i r,,ir)U [B^r') - e^i + ^.t] Ci(r')^.T(r') (28) 



Therefore 

fi(r, r', CO) = ^ J(r, r', co) + 4,(r) J^ ^^ f^t ^'^^ ' '^^^ ^;T(r)^H(r)Ci(r')V^.T(r') 

(29) 
where an expression for frequency dependent exchange interactions has the following form 

J(r, r', uj)=^-J^ ^'\~t\ r,,ir)B.c{r)^Hirmi{r')BUr')^,^{r') (30) 



l^u 



The later coincides with the exchange integrals |l|, y, |25[ if we neglect the u - dependence. 
Since B^^ ^ m we have J ~ m^ and the expression ()29|) vanishes in non- magnetic case, as 
it should be. Using the identity (j22j) one can show that 

f^(r, r', 0) = ^ J(r, r', 0) - B,,{r)S{r - r') (31) 

m[r) 

Note that for cj = we have exactly: 

and the static susceptibility x^^ (0) can be represented in the form 

£+- (0) = m (h-' - B-A (32) 



which is equivalent to the result of Ref. |l3 



n = h (^1 - B-^'hy' (33) 

for the renormalized exchange interaction if one define them in terms of inverse static sus- 
ceptibility liiLsl. 

As it was stressed above for a generic case of an itinerant electron magnet it is impossible 
to introduce the effective exchange integrals and one should to work with the generalized spin 
susceptibility. Any definition of the exchange integrals assume the adiabatic approximation 
somewhere. For the spin-wave spectrum which is determined by the pole of the transverse 
susceptibility it is natural to formulate "exchange concept" as neglecting of the LJ-dependence 
in Q. Then, in virtue of Eq. (j26p the magnon frequencies are just eigenstates of the operator 
Q (0) which exactly corresponds to the expression from the "old" MFT exchange interactions 
|l|, |2]. Note that for the long- wavelength limit q ^ this result turns out to be exact which 
proves the above statement about the stiffness constant D: in the framework of the local 



approximation it is rigorous. Corrections to D from a nonlocality of the exchange-correlation 
potential have been estimated recently J26| for Fe and Ni and turned out to be small. 

At the same time, if we are interested in the computations of the thermodynamic prop- 
erties such as the Curie temperature Tq the renormalized exchange integrals can really give 
more accurate results. One can introduce for the itinerant electron magnets a magnon-like 
operators 

6q = -L^-, 6i = -L5+ (34) 



where mo = 2S is the ground-state magnetization and write for the temperature dependence 
of the magnetization the Bloch-like expression 



oo 



m (T) = mo - V (6],&q) = mo + — V / du ' (35) 



q q 



(a similar approximation for the Hubbard model has been proposed in Ref. 1241 ). If we 
will use the classical-spin approximation usually exploited for first-principle estimations of 
the Curie temperature one should replace the Planck function in Eq. ()35p by its classical 
limit T /uj which immediately gives (taking into account the Kramers-Kronig relations) the 
following expression for the Curie temperature 

Tc m'o^ 

which is identical to the expression from Ref. |l3| in terms of the renormalized exchange 
interactions. Note, however, that the quantum character of the spin (which can be taken 
into account only beyond the LSDA) is probably very essential for proper description of high- 



temperature magnetism of transition metals [27[ which makes a problem of an improvement 
of classical estimations of Tc less important. 

In order to test different approximations to the exchange interactions we calculated spin- 



wave spectrum for iron and nickel using the LMTO-TB method 
representation was used and the calculation scheme was the following. 
The matrix of Green function in s,p,d-basis set is equal to 



28|. The orthogonal LMTO 



^..(k, Un) = [lUn + fi- Hl^Ai^)]-' (37) 



where /z is the chemical potential, uon are Matsubara frequencies and Hlda is the orthogonal 
LSDA Hamiltonian. We use the following approximation for the matrix Ixc'- 



V,, ^ ml,, = Hl^^iO) - Hi^^iO) 



(38) 



The matrix of the LDA-susceptibility has been calculated using the fast Fourier transform 
technique with k = (k, Un) 



Xr(g) = -E^T(A:)*Gi(A: + g) 

k 

The Fourier transforms of the "bare" exchange interactions J (0) (|3(J|) is defined as 



(39) 



^(q) = 4K.cXo+"(q,^ = 0)K, 



whereas for the "renormalized" exchange integrals 



la 



m 



one has 



J(0) — J(q) = -Ttl m {x^ {(\,uj = 0)) m 



-Trj 



m 



,-1 



{xt~ici,uj = 0))- -(x+-(q = 0,a; = 0)) 



-1 



m 



The magnon spectrum is determined via the exchange integrals as 



^(q) = ^ [^(0) - J(q)] 



(40) 



(41) 



(42) 



where M = Trim is the total magnetic moment. One can see from Fig.l and Fig. 2 that the 
LDA "bare" exchange parameter better describe the spin-wave spectrum in Fe and Ni, while 
thermodynamics (e.g. the Curie temperature - see Table. I) are more reasonable with the 
exact static LDA-exchange (which is the "RPA"-like expression). This confirms a general 
consideration presented above. 
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FIG. 1: The spin- wave spectrum for ferromagnetic iron in the bare exchange (Eg . (1401) ) and renor- 
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FIG. 2: The spin- wave spectrum for ferromagnetic nickel in the bare exchange fEa. (|4UI) ') and 
renormalized exchange (Ea.H41[)) scheme in comparison with experimental data (from Ref.j2 
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